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Abstract 

The current experimental and theoretical status of hadron and jet production at large transverse momentum 
in high-energy nucleus-nucleus collisions is summarised. The most important RHIC results are compared to 
theoretical parton energy loss predictions providing direct information on the (thermo)dynamical properties 
of hot and dense QCD matter. Prospects for the LHC are also outlined. 



1. Introduction 

The physics programme of high-energy nucleus-nucleus (AA) collisions is focused on the 
study of the fundamental theory of the strong interaction - Quantum Chromo Dynamics (QCD) 
- in extreme conditions of temperature, density and small parton momentum fraction (low-x) - 
see e.g. |1| for a recent review. By colliding two heavy nuclei at ultrarelativistic energies one 
expects to form a hot and dense deconfined medium whose collective (colour) dynamics can 
be studied experimentally. Lattice QCD calculations ^ predict a new form of matter at energy 
densities (well) above e„i, w 1 GeV/fm^ consisting of an extended volume of deconfined and 
chirally-symmetric (bare-mass) quarks and gluons: the Quark Gluon Plasma (QGP). 

One of the first proposed "smoking guns" of QGP formation was jet quenching i.e. the 
attenuation or disappearance of the spray of hadrons resulting from the fragmentation of a parton 
having suffered energy loss in the dense plasma formed in the reaction (Fig. [T] left). The en- 
ergy lost by a parton provides "tomographic" information of the matter properties (temperature 
T, interaction coupling a, thickness L): AE = f{E;T,a,L) |8|. The "scattering power" of the 
medium is often encoded in the transport coefficient which describes the average transverse mo- 
mentum squared transferred to the traversing parton per unit path-length: q = mj^/X ~ mjj p a 
(here mo is the medium Debye mass, p its density, and a the parton-matter interaction cross 
sectionQ The dominant mechanism of energy loss of a fast parton in a dense QCD plasma is 
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Fig. 1 . Right: "Jet quencliing" in a liead-on lieavy-ion collision: a fast parton traverses the dense plasma created (with 
transport coefficient q, gluon density dN^ /dy and temperature T), loses energy via "gluonstrahlung" and fragments into a 
(quenched) jet (3- Left: Neutral pion spectrum measured by PHENIX at ■Jsnn = 200 GeV in central AuAu (squares) (5], 
compared to the (T^A-scaled) spectrum in pp collisions (circles) [6] and to a NLO pQCD calculation (yellow band) (T)- 



of radiative nature ("gluonstrahlung"): the parton loses energy mainly by medium-induced mul- 
tiple gluon emission IIIOII 1I12I13I . Jet quenching in AA reactions is characterised by various 
observable consequences compared to the same "QCD vacuum" measurements in proton-proton 
{pp) collisions: (i) suppressed high-pj hadron spectrum (dNAA/dpr), (ii) unbalanced back-to- 
back high-pT dihadron azimuthal correlations {dNpair/dif), and (iii) modified energy -particle 
flow (softer hadron spectra, larger multiplicity, increased angular broadening, ...) within the final 
jets. A detailed review of these topics can be found in |4|, of which a summary is given in the 
following sections. 



2. High-pr single inclusive hadron production 

If a hard scattered parton suffers energy loss in a heavy-ion collision, the energy available for 
the hadrons issuing from its fragmentation will be reduced and their spectrum depleted compared 
to pp collisions. The standard method to quantify the medium effects on the yield of a large-pr 
particle produced at rapidity y in a AA reaction is given by the nuclear modification factor. 

RAA{pT,y;b) = n^^^it'^^^], , — 1 TaaI^) being the nuclear overlap function at Zj, (1) 
{TAA[b)) X d^Opp/dydpT 

which measures the deviation of AA at impact parameter b from an incoherent superposition of 
nucleon-nucleon collisions (Raa = !)■ From the measured suppression factor one can determine 
various medium properties such as its transport parameter q, via {AE) °^ ttj (q) L? Ul.lSJ . or its 
initial gluon density dN^ /dy, via AE ^ ajCRj-^'-^L (for an expanding plasma with original 

transverse area Aj_ = nRj^ w 150 fm^ and thickness L) lfT2l . We summarise the main high-/?7- 
hadroproduction results in pp and AA collisions, and confront them to jet quenching predictions. 
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(a) Magnitude of the suppression and medium properties.Figure[T] (right) shows the high-pj- 
Tt" spectrum measured at ^/s^ = 200 GeV in central AuAu [S"] compared to the pp (6\ and NLO 
pQCD tZJ spectra scaled by Taa- The AuAu data are suppressed by a factor of 4 - 5 with re- 
spect to the pp results. The corresponding Raa{pt), Eq. 0]), is shown in Fig. |2] (left). Above 
Pt ~ 5 GeV/c, Tt" f T4l . T] iQ?!, and charged hadrons I116I17II show all a common factor of ^5 
suppression relative to the Raa = 1 expectation that holds for hard probes, such as direct pho- 
tons 0181191 . which do no interact with the medium. The AuAu high-/??' suppression can be well 



"10 - 




Au+Au - 200 GeV (central collisions): 
■ □ Direct T, T* [PHENIX] 

* Inclusive h' [STAR] 
A It" [PHENIX] 

• ^ [PHENIX] 
GLV energy loss (dNf /dy = 1 400) 



10-^ 



^art Scaling 



•■^-TP-ij.^ — — 



8 10 12 14 16 18 20 

(GeV/c) 




10" 



High p in centrai AA: 

PbPb \/s^= 17.3 GeV [WA9B) 
■ AuAu \|s^ = 200 GeV (PHENiX) 
GLV energy ioss (dN /dy = 400) 
GLV energy ioss (dN /dy = 1400) 



RHIC 



GLV: dN(/dy = 2000 - 4000 

LHC 




PQM:<q> = 30-80GeV'^/lm 



10 



10 p^(GeV/c) 



Fig. 2. Left: Raa{pt) in central AuAu at 200 GeV for It" |5 |, r] |"T5l, charged hadrons fT6l, and direct y 1 18 19 | compared 
to the GLV model (dN^ /dy = 1400, yellow curve) |20|. Right: Raa{pt) for 7t"'s at SPS 121.22.1 and RHIC [5J compared 
to GLV calculations (dN^ jdy = 400, 1400) and to predictions for central PbPb at = 5.5 TeV (yellow bands) (23]: 

GLV (dNS /dy = 2000 - 4000) and PQM ({g) « 30 - 80 GeV^/fm). 



reproduced by parton energy loss calculations in a very dense system with initial gluon rapidity 
densities dN^ /dy w 1400 (Gyulassy-Levai-Vitev curve in Fig. 121 left) f20l, transport coefficients 
^ « 13 GeV^/fm [5.24J . or plasma temperatures T w 0.4 GeV [25 1. The consistency between the 
extracted q, dN^ /dy and T values in the various models has been studied e.g. in 041261 . Whereas 
the agreement between the fitted thermodynamical variables dN^ /dy and T is good, the values 
of the transport parameter q favoured by the data are 3 - 4 times larger than perturbative esti- 
mates |l9l. An accord between the obtained q and dN^ /dy can only be achieved assuming parton- 
medium cross-sections much larger than the <3gg « 1 .5 mb LO expectation. Such an observation 
lends support to the strongly-coupled nature of the QGP produced at RHIC 1.27 j . 



(b) Centre-of-mass energy dependence. As one increases the collision energy in nucleus- 
nucleus collisions, the produced plasma reaches higher energy and particle densities, the system 
stays longer in the QGP phase, and correspondingly the traversing partons are more quenched. 
Figure |2] (right) compiles the measured Raa{pt) for high-y^j- Ti" in central AA collisions in the 
range y/s///^ « 20 - 200 GeV compared to parton energy loss calculations that assume the forma- 
tion of a QGP with initial gluon densities in the range dN^ /dy k, 400 - 1400 II20I28II or, equiva- 
lently, averaged transport coefficients (^) w 3.5 - 13 GeV^/fm |24|. The theoretical predictions 
reproduce well the magnitude and shape of the experimental data. The SPS Raa (pt) Ii211 , though 
consistent with unity [22 1, is suppressed compared to the "Cronin enhancement" observed in pe- 
ripheral PbPb and pPb collisions at the same ^/s/^ [|29l . 
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(c) /?7--dependence of the suppression. At RHIC top energies, the hadron quenching factor 
remains basically constant from 5 GeV/c up to the highest transverse momenta measured so far, 
PT w 20 GeV/c (Fig. ^. Full calculations 1201241301311 including the combined effect of (i) 
energy loss kinematics constraints, (ii) steeply falling pr spectrum of the scattered partons, and 
(iii) ^(20%) pr-dependent (anti)shadowing differences between the proton and nuclear parton 
distribution functions (PDFs), result in an effectively flat Raa (pt) as found in the data. The much 
larger kinematical range opened at LHC energies 1231 will allow one to test the pj-dependence 
of parton energy loss over a much wider domain than at RHIC (yellow bands in Fig.|2l right). 

(d) Centrality (system-size) dependence. The volume of the produced plasma in a heavy-ion 
collision can be "dialed" by modifying the overlap area between the colliding nuclei either by 
selecting a given impact-parameter b - i.e. by choosing more central or peripheral reactions - or 
by colliding larger or smaller nuclei, e.g. Au (A = 197) versus Cu (A = 63). The relative energy 
loss depends on the effective mass number Aeff or, equivalently, on the number of participant 
nucleons in the collision A^pan, as: AE/E A^g o= 1241321 . The measured Raa{pt) in 
central CuCu at 22.4, 62.4, and 200GeV |331 is a factor of {AAu/Acufl^ « 2 lower than in 
central AuAu at the same energies. Yet, for a comparable Npan value, the suppression in AuAu 
and CuCu is very similar. Fitting the Npart dependence to Raa = (1 ^ '^^pan)"^^ yields a = 
0.56 ± 0.10 (5], consistent with parton energy loss calculations. 

(e) Path-length dependence.Experimentally, one can test the dependence of parton suppression 
on the plasma thickness (L) by exploiting the spatial asymmetry of the system produced in non- 
central nuclear collisions. Partons produced "in plane" ("out-of -plane") i.e. along the short (long) 
direction of the ellipsoid matter with eccentricity 8 will comparatively traverse a shorter (longer) 
thickness. PHENIX has measured the high-pj- neutral pion suppression as a function of the angle 
with respect to the reaction plane, Raa {pt , f34l . Each azimuthal angle (|) can be associated with 
an average medium path-length Lg via a Glauber model. The energy loss is found to satisfy the 
expected AE L dependence [12] above a "threshold" length of L w 2 fm. 



■ Pb+Pb ^ 7t"+X 0-77o central [WA98] 




x/Sn,; (GeV) Pt [GeV] 

Fig. 3. Left: Raa (pj = 4 GeV/c) for Tl" in central AA collisions as function of collision energy compared to non-Abelian 
(solid) and "non-QCD" (dotted) energy loss curves [35^36]. Right: Raa{pt) for decay electrons from D and B mesons in 
central AuAu at ^Js^n = 200 GeV I37I38I39I compared to a model of D and B meson dissociation in the plasma 1401 . 
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(f) Non-Abelian (colour factor) dependence. The amount of energy lost by a parton in a 
medium is proportional to its colour Casimir factor: Ca = 3 for gluons, Cp = 4-/3 for quarks. 
Asymptotically, the probability for a gluon to radiate another gluon is Ca/Cf = 91 A times larger 
than for a quark and thus g-jets are expected to be more quenched than ^-jets in a QGP. One 
can test such a genuine non-Abelian property of QCD energy loss by measuring hadron sup- 
pression at & fixed pt for increasing ^/s [35,361 . At large (small) x, the PDFs are dominated by 
valence-quarks ("wee" gluons) and consequently hadroproduction is dominated by quark (gluon) 
fragmentation. Figure|3](left) shows the Raa for4-GeV/c pions measured at SPS and RHIC com- 
pared to two parton energy loss curves l(36J . The lower (upper) curve shows the expected Raa 
assuming a normal (arbitrary) behaviour with AEg/AEg = 9/4 (AEg = AEg). The experimental 
high- Pt ti" data support!^ the expected colour-factor dependence of Raa ( ^/smi) lf35l . 



(g) Heavy-quark mass dependence. Due to the "dead cone" effect ll42ll . the radiative energy 
loss for a charm (bottom) quark is a factor l-(mg/mD) « 25% (75%) smaller than for a light- 
quark. Yet, RHIC measurements II37I38I39]| of high-/?/- electrons from the semi-leptonic decays 
of D- and B-mesons (Fig. [3] right) indicate the same suppression for light and heavy mesons: 
Raa{Q) ^ RAA{q,g) ~ 0.2, in contradiction with parton energy loss predictions 1431441 . Vari- 
ous explanations have been proposed to solve the 'heavy flavor puzzle' H. Among them is the 
observation that the hypothesis of vacuum hadronisation (after in-medium radiation) implicit in 
all formalisms may well not hold in the case of heavy quarks. The formation time of D- and 
B-mesons is of order T^^^.^^^ w 0.4 - 1 fm/c respectively and, thus, one needs to account for the 
energy loss of the heavy-quark as well as the possible dissociation of the heavy-quark meson 
inside the QGP |40|. The expected amount of suppression in that case is larger and consistent 
with the data (Fig. [5] right). 

3. High-/7r di-hadron correlations 

Jet-like correlations in heavy-ion collisions can be measured on a statistical basis by selecting 
a high- Pt trigger particle and measuring the azimuthal (A(|) = (|) — ^trig) and pseudorapidity (Ar| = 

^~^trig) distributions of its aiioc/afeii hadrons (/?^"™' < pj'^): C(A(|). Ar|) — dM)dKx\ • PP 
collisions, a dijet signal appears as two distinct back-to-back Gaussian-like peaks at A(|) w 0, 
Arj K, (near-side) and at A(|) w 7i (away-side). At variance with such a topology. Fig. |4] shows 
the increasingly distorted back-to-back azimuthal correlations in high-pj- triggered central AmAm 
events as one decreases the pT of the associated hadrons (right to left). Whereas, the AmAm and pp 
near-side peaks are similar for all pr's, the away-side peak is only present for the highest partner 
pr's but progressively disappears for less energetic partners II46I45I . Early STAR results BtI 
showed a monojet-like topology with a complete disappearance of the opposite-side peak for 
^a.«o6 ^2-4 GeV/c. The correlation strength over an azimuthal range A(|) between a trigger 
hadron /i, and a partner hadron ha in the opposite azimuthal direction can be constructed as a 
function of the momentum fraction zt = pj^""^ / p'j^ via a "pseudo-fragmentation function": 



The use of high-pr (anti)protons (mostly coming from gluon fragmentation) as an alternative test of the colour charge 
dependence of the quenching |4r| is, unfortunately, distorted by the presence of extra non-perturbative mechanisms of 
baryon production (see discussion in [4)). 
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Fig. 4. Comparison of the azimuthal di-hadron correlation dN p„ir / dA^d\\ for pp (open symbols) and central AuAu (his- 
tograms) at sJsnn = 200 GeV for pj'^ = 5-10 GeV/c and increasingly smaller (right to left) values of pf""^ 1451 . 



izr) - J dp, J dp, J dA<^ . (2) 

shown in Fig. |5] (top-left) compared to predictions of the HT model for various values of the 
£o parameter quantifying the amount of energy loss f^Sl. Similarly to Raa, the magnitude of 
the suppression of back-to-back jet-like two-particle correlations can be quantified with the ratio 
jcmay jymay j jyiway ^ joway fg^j^j jq decrease with increasing centrality, down to about 0.2 - 
0.3 for the most central events (Fig.|5] bottom-left) 1147 149 II . The right plot of Fig. |5] shows the 




best Eo ~ 1.9 GeV/fm value that fits the measured Raa{pt) and Iaa{zt) factors. Due to the 
irreducible presence of (unquenched) partons emitted from the surface of the plasma, the single- 
hadron quenching factor /?aa(/'7') is in general less sensitive to the value of £o than the dihadron 
modification ratio Iaa {zt ) ■ 
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4. Jet observables 



The measurements in collisions of fully reconstructed (di)jets or of jets tagged by an away- 
side photon or Z-boson allow one to investigate - in much more detail than with single- or double- 
hadron observables - the mechanisms of in-medium parton radiation as well as to characterise 
the matter properties through modified jet profiles 150 51] and fragmentation functions 11521 . 
Experimental reconstruction of jets in nuclear reactions is an involved three-steps exercise ||4l: 

(i) Clustering algorithm: The measured hadrons are clustered together, according to relative 
distances in momentum and/or space, following an infrared- and collinear-safe procedure 
which is also/flif enough to be run over events with very high multiplicities. The kj and 
SISCone algorithms implemented in the FastJet package |53| fulfill all such conditions. 

(ii) Background subtraction: Jets are produced on top of a large "underlying event" (UE) of 
hadrons coming from other (softer) parton-parton collisions in the same interaction. In 
central PbPb collisions at the LHC one expects Ej^ « 80 GeV (with large fluctuations) in 
a cone of radius R = y^rp+A(p = 0.4. Various UE subtraction techniques are available 
in combination with the kj ll54l . UAl-cone [55^| or iterative-cone |56| algorithms. 

(iii) Jet energy corrections: Data-driven methods are needed to experimentally control the jet 
energy-scale which is the single most important source of systematic uncertainties in the jet 
yield. The non-perturbative effects introduced by the UE and hadronisation can be gauged 
by comparing the sensitivity of the jet spectrum obtained with different Monte Carlo's ijU . 

STAR |[57| has a preliminary measurement of jets in AuAu at 200 GeV (Fig. |6] left) using a cone 
algorithm with R = 0.4, and estimating the UE background from the average energy in cones 
without seeds. Control of the jet energy corrections is still work in progress. Jet physics will 
definitely benefit from the highest energies (and therefore statistics) available at the LHC. The 
expected pj reach in PbPb at 5.5 TeV is as large as pj ~ 500 GeV/c (Fig.|6] right). 
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Fig. 6. Left: Preliminary STAR jet Ej .spectra in central AuAu (triangles) and pp (squares, scaled by Taa) at 200 GeV 1571 . 
Right: Jet spectra for various PhPb centralities expected at 5.5 TeV in CMS (fj^dt = 0.5 nb"') I58l . 

The y-jet (and Z-jet) channel provides a very clean means to determine medium-modified 
parton fragmentation functions (EEs) [ 59.60J . Since the prompt y is not affected by final-state 
interactions, its transverse energy (Ej) can be used as a proxy of the away-side parton energy 
(Ej' « E^) before any jet quenching. The EE, defined as the normalised distribution of hadron 
momenta 1 /NjeisdN/dz relative to that of the parent parton E^', can then be constructed using 
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■^y;, ~ Pt /Ej' or, similarly, ^ — \og{Ej / pj) = — log(z), for all particles with momentum pj 
associated with a jet. In a QCD medium, energy loss shifts parton energy from high-z to low-z 
hadrons 161], resulting in a higher "hump-back plateau" in the FFs at intermediate ^ w 3 - 4 
values (Fig. |7] left). Full simulation studies of the y-jet channel in central Pb-Pb (Fig. |7] right) 
indicate that medium modified FFs are measurable with small uncertainties in the ranges z < 0.7 
and 0.2 <£ < 6 [62J. 
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Fig. 7. Left: Single inclusive distribution of liadrons vs. ^ = \n{Ejg,jp) for a 17.5-GeV jet in e+e" collisions (TASSO 
data) compared to QCD radiation predictions in the vacuum (solid curve) and in-medium (dashed curve) 1611 . Right: FFs 
as a function of ^ for quenched partons obtained in CMS y-jet simulations for central Pb-Pb at 5.5 TeV (0.5 nb" ' ) 1 621 . 



5. Summary 

The analysis of jet structure modifications in heavy-ion collisions provides quantitative "tomo- 
graphic" information on the thermodynamic al and transport properties of the strongly interacting 
medium produced in the reactions. At RHIC energies (up to y/sj^ = 200 GeV), strong suppres- 
sion of the yields of high-pr single hadrons and of dihadron azimuthal correlations, have been 
observed in central AuAu collisions. Most of the properties of the observed suppression are in 
quantitative agreement with the predictions of parton energy loss models in a very dense QCD 
plasma. The confrontation of these models to the data permits to derive the initial gluon density 
dN^ /dy w 1400 and transport coefficient q— ff{10 GeV^/fm) of the produced medium at RHIC. 
At the upcoming LHC energies, the detailed analysis of jet spectra, jet shapes and the extrac- 
tion of medium-modified parton-to-hadron fragmentation functions promise to fully unravel the 
mechanisms of parton energy loss in QCD matter. The study of jet quenching phenomena proves 
an excellent tool to expand our knowledge of the dynamics of the strong interaction at extreme 
conditions of temperature and density. 



Acknowledgments 

Special thanks to Itzhak Tserruya and the organisers of PANIC 08 for their kind invitation. Sup- 
port by the 6th EU Framework Programme contract MEIF-CT-2005-025073 is acknowledged. 



8 



References 



[I] D. d'Enterria, J. Phys. G 34, S53 (2007) 

[2] M. Cheng et al, Phys. Rev. D 74, 054507 (2006); Y. Aoki et al, Phys. Lett. B 643, 46 (2006) 

[3] J. D. Bjorken, FERMILAB-PUB-82-059-THY (1982) 

[4] D. d'Enterria, arXiv:0902.2011 [nucl-ex] 

[5] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 101, 232301 (2008) 

[6] A. Adare et al. [PHENIX Collaboration], Phys. Rev. D 76, 05 1 106 (2007) 

[7] F. Aversa, P Chiappetta, M. Greco and J. P Guillet, Nucl. Phys. B 327, 105 (1989); B. Jager, A. Schafer, 

M. Stratmann and W. Vogelsang, Phys. Rev. D 67, 054005 (2003); W. Vogelsang (private communication) 

[8] S. Peigne and A. V. Smilga, ai-Xiv:08 10.5702 [hep-ph] 

[9] R. Baier and D. Schiff, JHEP 0609, 059 (2006) 

[10] M. Gyulassy, M. Plumer, Phys. Lett. B243, 432 (1990); X.N. Wang, M. Gyulassy, Phys. Rev. Lett. 68, 1480 (1992) 

[I I] R. Baier, Y.L. Dokshitzer, A.H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B484, 265 (1997); R. Baier, D. Schiff, 
B.G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50, 37 (2000) 

[12] M. Gyulassy, P Levai and 1. Vitev, Phys. Rev. Lett. 85, 5535 (2000); Nucl. Phys. B594, 371 (2001) 

[13] U. A. Wiedemann, Nucl. Phys. B 588, 303 (2000) 

[14] S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91, 072301 (2003) 

[15] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 96, 202301 (2006) 

[16] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 91, 172302 (2003) 

[17] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C69, 034910 (2004) 

[18] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 94, 232301 (2005) 

[19] A. Adare e( a/. [PHENIX Collaboration], arXiv:0804.4168 [nucl-ex] 

[20] 1. Vitev and M. Gyulassy, Phys. Rev. Lett. 89, 252301 (2002); 1. Vitev, J. Phys. G 30, S791 (2004) 

[21] M. M. Aggarwal et al. [WA98 Collaboration], Eur. Phys. J. C 23, 225 (2002) 

[22] D. d'Enterria, Phys. Lett. B 596, 32 (2004) 

[23] N. Armesto et al, J. Phys. G 35, 054001 (2008) 

[24] A. Dainese, C. Loizides and G. Paic, Eur. Phys. J. C 38, 461 (2005) 

[25] S. Turbide, C. Gale, S. Jeon and G. D. Moore, Phys. Rev. C 72, 014906 (2005) 

[26] S. A. Bass etal, arXiv:0808.0908 [nucl-th] 

[27] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 (2005) 

[28] 1. Vitev, Phys. Lett. B 606, 303 (2005) 

[29] M. M. Aggarwal et al. [WA98 Collaboration], Phys. Rev. Lett. 100, 242301 (2008) 

[30] S. Jeon and G. D. Moore, Phys. Rev C 71, 034901 (2005) 

[31] K. Eskola, H. Honkanen, C. Salgado and U. Wiedemann, Nucl. Phys. A 747, 5 1 1 (2005) 

[32] 1. Vitev, Phys. Lett. B 639, 38 (2006) 

[33] A. Adai-e et al. [PHENIX Collaboration], Phys. Rev. Lett. 101, 162301 (2008) 

[34] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 76, 034904 (2007) 

[35] D. d'Enterria, Eur. Phys. J. C43, 295 (2005) 

[36] Q. Wang and X.N. Wang, Phys. Rev. C71, 014903 (2005) 

[37] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 96, 032301 (2006) 

[38] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 172301 (2007) 

[39] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 98, 192301 (2007) 

[40] A. Adil and 1. Vitev, Phys. Lett. B 649, 139 (2007) 

[41] B. Mohanty [STAR Collaboration], J. Phys. G 35, 104006 (2008) 

[42] Y. L. Dokshitzer and D. E. Khai'zeev, Phys. Lett. B 519, 199 (2001) 

[43] M. Djordjevic, M. Gyulassy and S. Wicks, Phys. Rev. Lett. 94, 112301 (2005) 

[44] N. Armesto, M. Cacciari, A. Dainese, C. A. Salgado and U. A. Wiedemann, Phys. Lett. B 637, 362 (2006) 

[45] A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 78, 014901 (2008) 

[46] A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 77, 01 1901 (2008) 

[47] C. Adler et al. [STAR Collaboration], Phys. Rev Lett. 90, 082302 (2003) 

[48] H. Zhang, J. F. Owens, E. Wang and X. N. Wang, Phys. Rev. Lett. 98, 212301 (2007) 

[49] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 97, 162301 (2006) 

[50] C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 93, 042301 (2004) 



9 



[51] I. Vitev, S. Wicks and B. W. Zhang. l^Civ:08l0.2807i rheD-ph1 
[52] F. Ai-leo, arXiv:0810.H93 [hep-ph] 
[53] M. Cacciari and G. P. Salam, Phys. Lett. B 641, 57 (2006) 
[54] M. Cacciari and G. P Salam, Phys. Lett. B 659, 1 19 (2008) 
[55] S. L. Blyth et ai, J. Phys. G 34, 271 (2007) 

[56] O. Kodolova, L Vardanian, A. Nikitenko and A. Oulianov, Eur. Phys. J. C 50, 117 (2007) 

[57] J. Putschke [STAR Collaboration], arXiv:0809.1419 , S. Salur [STAR Collaboration], arXiv:0810.0500l lnucl- 

[58] D. d'Enterria (ed.) et al. [CMS Collaboration], J. Phys. G 34, 2307 (2007) 

[59] X. N. Wang, Z. Huang and 1. Sarcevic, Phys. Rev. Lett. 77, 231 (1996) 

[60] F. Arleo, P Aurenche, Z. Belghobsi and J. P Guillet, JHEP 0411, 009 (2004) 

[61] N. Borghini and U. A. Wiedemann, arXiv:hep-ph/0506218 

[62] C. Loizides [CMS Collaboration], J. Phys. G 35, 104166 (2008) 



10 



